Exploring the path to polio eradication: insights from consecutive seroprevalence surveys among Pakistani children.

Introduction: After trivalent oral poliovirus vaccine (tOPV) cessation, Pakistan has maintained immunity to type 2 poliovirus by administering inactivated polio vaccine (IPV) in routine immunization, alongside monovalent OPV type 2 (mOPV2) and IPV in supplementary immunization activities (SIAs). This study assesses the change in poliovirus type 2 immunity after tOPV withdrawal and due to SIAs with mOPV2 and IPV among children aged 6-11 months. Methods: Three cross-sectional sequential serological surveys were conducted in 12 polio high-risk areas of Pakistan. 25 clusters from each geographical stratum were selected utilizing probability proportional to size. Results: Seroprevalence of type 2 poliovirus was 49%, with significant variation observed among surveyed areas; <30% in Pishin, >80% in Killa Abdullah, Mardan & Swabi, and Rawalpindi. SIAs with IPV improved immunity from 38 to 57% in Karachi and 60 to 88% in Khyber. SIAs with IPV following mOPV2 improved immunity from 62 to 65% in Killa Abdullah, and combined mOPV2 and IPV SIAs in Pishin improved immunity from 28 to 89%. Results also reflected that immunity rates for serotypes 1 and 3 were consistently above 90% during all three phases and across all geographical areas. Conclusion: The study findings highlight the importance of implementing effective vaccination strategies to prevent the re-emergence of poliovirus. Moreover, the results provide crucial information for policymakers working toward achieving global polio eradication.

Modelling the spread of serotype-2 vaccine derived-poliovirus outbreak in Pakistan and Afghanistan to inform outbreak control strategies in the context of the COVID-19 pandemic.

BACKGROUND: Since July 2019, Pakistan and Afghanistan have been facing an outbreak of serotype-2 circulating vaccine derived poliovirus (cVDPV2) in addition to continued transmission of serotype-1 wild poliovirus (WPV1) and SARS-CoV-2 in 2020. Understanding the risks of cVDPV2 transmission due to pause of global vaccination efforts and the impact of potential vaccination response strategies in the current context of COVID-19 mitigation measures is critical. METHODS: We developed a stochastic, geographically structured mathematical model of cVDPV2 transmission which captures both mucosal and humoral immunity separately and allows for reversion of serotype-2 oral polio vaccine (OPV2) virus to cVDPV2 following vaccine administration. The model includes geographic heterogeneities in vaccination coverage, population immunity and population movement. The model was fitted to historic cVDPV2 cases in Pakistan and Afghanistan between January 2010-April 2016 and July 2019-March 2020 using iterated particle filtering. The model was used to simulate spread of cVDPV2 infection from July 2019 to explore impact of various proposed vaccination responses on stopping transmission and risk of spread of reverted Sabin-2 under varying assumptions of impacts from COVID-19 lockdown measures on movement patterns as well as declines in vaccination coverage. RESULTS: Simulated monthly incidence of cVDPV2 from the best-fit model demonstrated general spatio-temporal alignment with observed cVDPV2 cases. The model predicted substantial spread of cVDPV2 infection, with widespread transmission through 2020 in the absence of any vaccination activities. Vaccination responses were predicted to substantially reduce transmission and case burden, with a greater impact from earlier responses and those with larger geographic scope. While the greatest risk of seeding reverted Sabin-2 was predicted in areas targeted with OPV2, subsequent spread was greatest in areas with no or delayed response. The proposed vaccination strategy demonstrated ability to stop the cVDPV2 outbreak (with low risk of reverted Sabin-2 spread) by February 2021. CONCLUSION: Outbreak response vaccination campaigns against cVDPV2 will be challenging throughout the COVID-19 pandemic but must be implemented urgently when feasible to stop transmission of cVDPV2.

Quantifying movement patterns and vaccination status of high risk mobile populations in Pakistan and Afghanistan to inform poliovirus risk and vaccination strategy.

BACKGROUND: Stopping serotype 1 wild poliovirus transmission in Pakistan and Afghanistan requires ensuring all children <5 years of age are repeatedly vaccinated, including the large proportion living in mobile groups. Vaccinating children living in high-risk mobile populations (HRMPs) remains a priority for the polio programme. METHODS: In 2017-2018, group-level censuses were conducted in 43 districts of Pakistan, gathering information for all HRMP children <5 years of age residing in settlements. Demographic and mobility information was collected, including HRMP type, ethnicity, language, mode of transportation and movement patterns. Vaccination status was recorded for the most recent polio campaign. Proportion of HRMP children by demographic factors and mode of transportation was determined and the magnitude of movement was quantified based on the origin, previous and next locations. Magnitude of cross-border movement with Afghanistan was evaluated, as was primary crossing point. Vaccination status was evaluated for each district by demographic and mode of transportation information. RESULTS: In total, 188,130 HRMP children <5 years of age were assessed. The predominant HRMP type, ethnic group, language and mode of transport was Afghan refugees (27%), Pashtun (69%), Pashto (69%) and bus (52%). Overall, 84% of children originated outside of their current district, including 29% from Afghanistan. Previous and next locations, were reported outside of current location by 34% and 77% of children. Afghanistan was previous and next location for 5% and 11% of children, with 5.5% and 3% of children crossing the Afghanistan border in the past 6-months and next 3-months. Primary crossing route was Torkham (79%). Overall vaccination coverage was 98% (IQR: 96%-99%) and consistently >90% across HRMP type, ethnic group, language and mobility means. CONCLUSION: Large numbers of HRMPs were found across Pakistan, with substantial links throughout the country and with Afghanistan. While vaccination coverage of HRMPs was high, ensuring these populations are consistently vaccinated remains a priority.

Routine immunization coverage in Pakistan: a survey of children under 1 year of age in community-based vaccination areas.

Pakistan is one of two countries in which poliovirus remains endemic. Considering the high number of children born every year, reaching and vaccinating new birth cohorts by improving routine immunization coverage in children <1 year of age is crucial to halting virus transmission. In 2015, a community-based vaccination (CBV) strategy, using local community members to enhance vaccine acceptance and improve routine immunization service delivery, was introduced in areas of Pakistan that have never interrupted poliovirus transmission. In order to assess progress towards improving routine immunization, we performed house-to-house immunization surveys across ten CBV areas in 2017 and 2018. In each household, we determined age-appropriate routine antigen coverage for children <1 year of age based on vaccination card and caregiver recall. We surveyed 5,499 and 5,264 children in 2017 and 2018, respectively. Overall, coverage of inactivated poliovirus vaccine (IPV) at 14 weeks of age was 32% in 2017 and 39% in 2018 based on vaccination card and recall. Across the surveyed areas, coverage ranged from 7% in Killa Abdullah to 61% in Peshawar in 2018. Oral poliovirus vaccination coverage decreased with successive vaccination visits, ranging from 66% for the birth dose to 42% for the 14-week dose in 2018. No area reached the target of 80% coverage for any routine antigen. Our findings highlight the need for concerted efforts to improve routine immunization coverage in these critical areas of wild poliovirus transmission.

Needle-free injectors for mass administration of fractional dose inactivated poliovirus vaccine in Karachi, Pakistan: A survey of caregiver and vaccinator acceptability.

The first large-scale vaccination campaign using needle-free jet injectors to administer fractional doses of inactivated poliovirus vaccine (fIPV) was conducted in Karachi, Pakistan, in February 2019. Data on acceptability of jet injectors were collected from 610 vaccinators and 4898 caregivers during the first four days of the campaign. Of those with prior needle and syringe experience, both vaccinators and caregivers expressed a strong preference for jet injectors (578/592 [97.6%] and 4792/4813 [99.6%], respectively), citing ease of use, appearance, and child's response to vaccination. Among caregivers, 4638 (94.7%) stated they would be more likely to bring their child for vaccination in a future campaign that used jet injectors. Mean vaccine coverage among towns administering fIPV was 98.7% - an increase by 18.4% over the preceding campaign involving full-dose IPV. Our findings demonstrate the strong acceptability of fIPV jet injectors and highlight the potential value of this method in future mass campaigns.

Assessing the sensitivity of the polio environmental surveillance system.

BACKGROUND: The polio environmental surveillance (ES) system has been an incredible tool for advancing polio eradication efforts because of its ability to highlight the spatial and temporal extent of poliovirus circulation. While ES often outperforms, or is more sensitive than AFP surveillance, the sensitivity of the ES system has not been well characterized. Fundamental uncertainty of ES site sensitivity makes it difficult to interpret results from ES, particularly negative results. METHODS AND FINDINGS: To study ES sensitivity, we used data from Afghanistan and Pakistan to examine the probability that each ES site detected the Sabin 1, 2, or 3 components of the oral polio vaccine (OPV) as a function of virus prevalence within the same district (estimated from AFP data). Accounting for virus prevalence is essential for estimating site sensitivity because Sabin detection rates should vary with prevalence-high immediately after supplemental immunization activities (SIAs), but low in subsequent months. We found that most ES sites in Pakistan and Afghanistan are highly sensitive for detecting poliovirus relative to AFP surveillance in the same districts. For example, even when Sabin poliovirus is at low prevalence of ~0.5-3% in AFP surveillance, most ES sites have ~34-50% probability of detecting Sabin. However, there was considerable variation in ES site sensitivity and we flagged several sites for re-evaluation based on low sensitivity rankings and low wild polio virus detection rates. In these areas, adding new sites or modifying collection methods in current sites could improve sensitivity of environmental surveillance. CONCLUSIONS: Relating ES detections to virus prevalence significantly improved our ability to evaluate site sensitivity compared to evaluations based solely on ES detection rates. To extend our approach to new sites and regions, we provide a preliminary framework for relating ES and AFP detection rates, and descriptions of how detection rates might relate to SIAs and natural seasonality.

Effect of Inactivated Poliovirus Vaccine Campaigns, Pakistan, 2014-2017.

Pakistan began using inactivated poliovirus vaccine alongside oral vaccine in mass campaigns to accelerate eradication of wild-type poliovirus in 2014. Using case-based and environmental surveillance data for January 2014-October 2017, we found that these campaigns reduced wild-type poliovirus detection more than campaigns that used only oral vaccine.

Population sensitivity of acute flaccid paralysis and environmental surveillance for serotype 1 poliovirus in Pakistan: an observational study.

BACKGROUND: To support poliomyelitis eradication in Pakistan, environmental surveillance (ES) of wastewater has been expanded alongside surveillance for acute flaccid paralysis (AFP). ES is a relatively new method of surveillance, and the population sensitivity of detecting poliovirus within endemic settings requires estimation. METHODS: Data for wild serotype 1 poliovirus from AFP and ES from January 2011 to September 2015 from 14 districts in Pakistan were analysed using a multi-state model framework. This framework was used to estimate the sensitivity of poliovirus detection from each surveillance source and parameters such as the duration of infection within a community. RESULTS: The location and timing of poliomyelitis cases showed spatial and temporal variability. The sensitivity of AFP surveillance to detect serotype 1 poliovirus infection in a district and its neighbours per month was on average 30.0% (95% CI 24.8-35.8) and increased with the incidence of poliomyelitis cases. The average population sensitivity of a single environmental sample was 59.4% (95% CI 55.4-63.0), with significant variation in site-specific estimates (median varied from 33.3-79.2%). The combined population sensitivity of environmental and AFP surveillance in a given month was on average 98.1% (95% CI 97.2-98.7), assuming four samples per month for each site. CONCLUSIONS: ES can be a highly sensitive supplement to AFP surveillance in areas with converging sewage systems. As ES for poliovirus is expanded, it will be important to identify factors associated with variation in site sensitivity, leading to improved site selection and surveillance system performance.

Routine immunization in Pakistan: comparison of multiple data sources and identification of factors associated with vaccination.

Background: Within Pakistan, estimates of vaccination coverage with the pentavalent vaccine, oral polio vaccine (OPV) and measles vaccine (MV) in 2011 were reported to be 74%, 75% and 53%, respectively. These national estimates may mask regional variation. The reasons for this variation have not been explored. Methods: Data from the Multiple Indicator Cluster Surveys (MICS) for Balochistan and Punjab (2010-2011) are analysed to examine factors associated with receiving three or more doses of the pentavalent vaccine and one or more MVs using regression modelling. Pentavalent and OPV estimates from the MICS were compared to vaccine dose histories from surveillance for acute flaccid paralysis (AFP; poliomyelitis) to ascertain agreement. Results: Adjusted coverage of children 12-23 months of age were estimated to be 16.0%, 75.5% and 34.2% in Balochistan and 58.0%, 87.7% and 72.6% in Punjab for the pentavalent vaccine, OPV and MV, respectively. Maternal education, healthcare utilization and wealth were associated with receiving the pentavalent vaccine and the MV. There was a strong correlation of district estimates of vaccination coverage between AFP and MICS data, but AFP estimates of pentavalent coverage in Punjab were biased toward higher values. Conclusions: National estimates mask variation and estimates from individual surveys should be considered alongside other estimates. The development of strategies targeted towards poorly educated parents within low-wealth quintiles that may not typically access healthcare could improve vaccination rates.

Fractional-Dose Inactivated Poliovirus Vaccine Campaign - Sindh Province, Pakistan, 2016.

Following the declaration of eradication of wild poliovirus (WPV) type 2 in September 2015, trivalent oral poliovirus vaccine (tOPV) was withdrawn globally to reduce the risk for type 2 vaccine-derived poliovirus (VDPV2) transmission; all countries implemented a synchronized switch to bivalent OPV (type 1 and 3) in April 2016 (1,2). Any isolation of VDPV2 after the switch is to be treated as a potential public health emergency and might indicate the need for supplementary immunization activities (3,4). On August 9, 2016, VDPV2 was isolated from a sewage sample taken from an environmental surveillance site in Hyderabad, Sindh province, Pakistan. Possible vaccination activities in response to VDPV2 isolation include the use of injectable inactivated polio vaccine (IPV), which poses no risk for vaccine-derived poliovirus transmission. Fractional-dose, intradermal IPV (fIPV), one fifth of the standard intramuscular dose, has been developed to more efficiently manage limited IPV supplies. fIPV has been shown in some studies to be noninferior to full-dose IPV (5,6) and was used successfully in response to a similar detection of a single VDPV2 isolate from sewage in India (7). Injectable fIPV was used for response activities in Hyderabad and three neighboring districts. This report describes the findings of an assessment of preparatory activities and subsequent implementation of the fIPV campaign. Despite achieving high coverage (>80%), several operational challenges were noted. The lessons learned from this campaign could help to guide the planning and implementation of future fIPV vaccination activities.

Spatial model for risk prediction and sub-national prioritization to aid poliovirus eradication in Pakistan.

BACKGROUND: Pakistan is one of only three countries where poliovirus circulation remains endemic. For the Pakistan Polio Eradication Program, identifying high risk districts is essential to target interventions and allocate limited resources. METHODS: Using a hierarchical Bayesian framework we developed a spatial Poisson hurdle model to jointly model the probability of one or more paralytic polio cases, and the number of cases that would be detected in the event of an outbreak. Rates of underimmunization, routine immunization, and population immunity, as well as seasonality and a history of cases were used to project future risk of cases. RESULTS: The expected number of cases in each district in a 6-month period was predicted using indicators from the previous 6-months and the estimated coefficients from the model. The model achieves an average of 90% predictive accuracy as measured by area under the receiver operating characteristic (ROC) curve, for the past 3 years of cases. CONCLUSIONS: The risk of poliovirus has decreased dramatically in many of the key reservoir areas in Pakistan. The results of this model have been used to prioritize sub-national areas in Pakistan to receive additional immunization activities, additional monitoring, or other special interventions.

Risk factors and short-term projections for serotype-1 poliomyelitis incidence in Pakistan: A spatiotemporal analysis.

BACKGROUND: Pakistan currently provides a substantial challenge to global polio eradication, having contributed to 73% of reported poliomyelitis in 2015 and 54% in 2016. A better understanding of the risk factors and movement patterns that contribute to poliovirus transmission across Pakistan would support evidence-based planning for mass vaccination campaigns. METHODS AND FINDINGS: We fit mixed-effects logistic regression models to routine surveillance data recording the presence of poliomyelitis associated with wild-type 1 poliovirus in districts of Pakistan over 6-month intervals between 2010 to 2016. To accurately capture the force of infection (FOI) between districts, we compared 6 models of population movement (adjacency, gravity, radiation, radiation based on population density, radiation based on travel times, and mobile-phone based). We used the best-fitting model (based on the Akaike Information Criterion [AIC]) to produce 6-month forecasts of poliomyelitis incidence. The odds of observing poliomyelitis decreased with improved routine or supplementary (campaign) immunisation coverage (multivariable odds ratio [OR] = 0.75, 95% confidence interval [CI] 0.67-0.84; and OR = 0.75, 95% CI 0.66-0.85, respectively, for each 10% increase in coverage) and increased with a higher rate of reporting non-polio acute flaccid paralysis (AFP) (OR = 1.13, 95% CI 1.02-1.26 for a 1-unit increase in non-polio AFP per 100,000 persons aged <15 years). Estimated movement of poliovirus-infected individuals was associated with the incidence of poliomyelitis, with the radiation model of movement providing the best fit to the data. Six-month forecasts of poliomyelitis incidence by district for 2013-2016 showed good predictive ability (area under the curve range: 0.76-0.98). However, although the best-fitting movement model (radiation) was a significant determinant of poliomyelitis incidence, it did not improve the predictive ability of the multivariable model. Overall, in Pakistan the risk of polio cases was predicted to reduce between July-December 2016 and January-June 2017. The accuracy of the model may be limited by the small number of AFP cases in some districts. CONCLUSIONS: Spatiotemporal variation in immunization performance and population movement patterns are important determinants of historical poliomyelitis incidence in Pakistan; however, movement dynamics were less influential in predicting future cases, at a time when the polio map is shrinking. Results from the regression models we present are being used to help plan vaccination campaigns and transit vaccination strategies in Pakistan.

Population Immunity against Serotype-2 Poliomyelitis Leading up to the Global Withdrawal of the Oral Poliovirus Vaccine: Spatio-temporal Modelling of Surveillance Data.

BACKGROUND: Global withdrawal of serotype-2 oral poliovirus vaccine (OPV2) took place in April 2016. This marked a milestone in global polio eradication and was a public health intervention of unprecedented scale, affecting 155 countries. Achieving high levels of serotype-2 population immunity before OPV2 withdrawal was critical to avoid subsequent outbreaks of serotype-2 vaccine-derived polioviruses (VDPV2s). METHODS AND FINDINGS: In August 2015, we estimated vaccine-induced population immunity against serotype-2 poliomyelitis for 1 January 2004-30 June 2015 and produced forecasts for April 2016 by district in Nigeria and Pakistan. Population immunity was estimated from the vaccination histories of children <36 mo old identified with non-polio acute flaccid paralysis (AFP) reported through polio surveillance, information on immunisation activities with different oral poliovirus vaccine (OPV) formulations, and serotype-specific estimates of the efficacy of these OPVs against poliomyelitis. District immunity estimates were spatio-temporally smoothed using a Bayesian hierarchical framework. Coverage estimates for immunisation activities were also obtained, allowing for heterogeneity within and among districts. Forward projections of immunity, based on these estimates and planned immunisation activities, were produced through to April 2016 using a cohort model. Estimated population immunity was negatively correlated with the probability of VDPV2 poliomyelitis being reported in a district. In Nigeria and Pakistan, declines in immunity during 2008-2009 and 2012-2013, respectively, were associated with outbreaks of VDPV2. Immunity has since improved in both countries as a result of increased use of trivalent OPV, and projections generally indicated sustained or improved immunity in April 2016, such that the majority of districts (99% [95% uncertainty interval 97%-100%] in Nigeria and 84% [95% uncertainty interval 77%-91%] in Pakistan) had >70% population immunity among children <36 mo old. Districts with lower immunity were clustered in northeastern Nigeria and northwestern Pakistan. The accuracy of immunity estimates was limited by the small numbers of non-polio AFP cases in some districts, which was reflected by large uncertainty intervals. Forecasted improvements in immunity for April 2016 were robust to the uncertainty in estimates of baseline immunity (January-June 2015), vaccine coverage, and vaccine efficacy. CONCLUSIONS: Immunity against serotype-2 poliomyelitis was forecasted to improve in April 2016 compared to the first half of 2015 in Nigeria and Pakistan. These analyses informed the endorsement of OPV2 withdrawal in April 2016 by the WHO Strategic Advisory Group of Experts on Immunization.